Nematic braids: topological invariants and rewiring of disclinations 
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The conventional topological description given by the fundamental group of nematic order pa- 
rameter does not adequately explain the entangled defect line structures that have been observed 
in nematic colloids. We introduce a new topological invariant, the self-linking number, that enables 
a complete classification of entangled defect line structures in general nematics, even without par- 
ticles, and demonstrate our formalism using colloidal dimers, for which entangled structures have 
been previously observed. We also unveil a simple rewiring scheme for the orthogonal crossing of 
two —1/2 disclinations, based on a tetrahedral rotation of two relevant disclination segments, that 
allows us to predict possible nematic braids and calculate their self-linking numbers. 



PACS numbers: 61.30.Dk,61.30.Jf,82.70.Dd 

Depending on the temperature and their molecular 
properties, nematogenic media can form isotropic, ne- 
matic, chiral nematic or even blue phases, each charac- 
terized by a specific orientational ordering of the con- 
stituent molecules. Their anisotropic nature allows the 
formation of disclinations that can be stabilized by geo- 
metric or intrinsic constraints. Recently, a lot of progress 
has been made on the stabilization and manipulation of 
disclinations by using dispersions of colloidal particles [Ti- 
ll] or confinement to porous networks [5] . In nematic dis- 
persions, the anisotropic inter-particle interactions medi- 
ated by elastic deformations and defects lead to diverse 
colloidal structures that promote self-assembly and of- 
fer great potential for photonics and plasmonics [6 . Ne- 
matic braids are disclination networks where defect loops 
are not localized around one particle, but instead entan- 
gle clusters of particles [3 [8]. The existence of entan- 
gled structures was first proposed based on the results 
of numerical simulations [9j HO] . They have since been 
observed experimentally and their stability has been ex- 
tensively analyzed [7J [11] . These structures are not suffi- 
ciently well described by the theory developed for simple 
nematic defects [12] and a complete theoretical under- 
standing is still lacking. Nematic braids may also include 
knots and links [13] , otherwise seen in the physics of poly- 
mers [14 , DNA [T5HT7] and knotted light [18 . Easy ex- 
perimental observation of nematic disclination networks, 
and their rewiring, knotting and linking by laser tweezers 
[3 5 places nematic braids as a primary template for 
the study of nontrivial topology in physical systems. 

Nematic braids stabilized by homeotropic particles 
consist of closed —1/2 disclination loops. To fully de- 
scribe a single disclination loop, we generalize the math- 
ematical notion of a loop by introducing the self-linking 
number, which counts how many times the cross section 
of the disclination turns during a complete loop. This 
invariant applies to elastic loops, DNA loops [17 and 
other fields, but in the case of a —1/2 nematic disclina- 
tion, due to its intrinsic three-fold symmetry, it assumes 
specific fractional values, similar to flux discretization in 
the fractional quantum Hall effect [T9] . 




FIG. 1. Rewiring sites of different dimer structures, (a-c) 
Theta structure and two chiral omega structures (simulations 
by M. Ravnik [7]). Rewiring sites are marked and paired with 
corresponding idealized structures, (d) The three conforma- 
tions of a disclination crossing with tetrahedral symmetry. 
Rewiring is performed by rotating around a C3 symmetry 
axis of the tetrahedron. 

For the investigation of the self-linking number, we 
chose a colloidal dimer consisting of two spherical par- 
ticles with strong homeotropic anchoring, confined to a 
homogeneous planar nematic cell [7]. Depending on the 
particle size, confinement and initial conditions, the par- 
ticles can interact by arranging themselves into dipolar 
or quadrupolar structures [2., or they can be bound by 
— 1/2 disclination loops shared between both particles 
[TT] , The homogeneous director field environment ener- 
getically disfavors linking and knotting, which reveals the 
more basic rewiring properties of entangled states. 

We demonstrate that differences between dimer struc- 
tures are localized to tetrahedral regions around cross- 
ings of disclinations (Fig. [I]), from which we derive rules 
for calculating the self-linking number and classifing all 
dimer structures. The rewiring rules apply to structures 
involving —1/2 disclinations in any confinement and can 
be used to predict and design nematic braids consisting 



of complex linked and knotted loops. 

We start by examining the similarities of two dimer 
disclination configurations [TT]. The "entangled hyper- 
bolic defect structure" (referred to as the theta struc- 
ture from here on) is the only dimer structure with space 
inversion symmetry and consists of two perpendicular 
loops, one encircling both particles and the other placed 
symmetrically between them (Fig. [IJi). The "omega 
structure", on the other hand, consists of a single loop 
wrapped around both particles and has two chiral isomers 
(Fig.[^3,c). All three structures have similar director field 
and line geometry at the far ends of the colloidal parti- 
cles and along the vertical axis (Fig. 5 in [H]). They only 
differ in the way the left arc, right arc and central loop 
meet between the particles (see the encircled areas on 
Fig.^c). The conversion of one structure into another 
is achieved by rewiring the crossing, while leaving the 
remote field intact, which requires a cutting of disclina- 
tions. The resulting four endpoints define a tetrahedron 
that encloses the rewiring site. Inside the tetrahedron, 
two perpendicular disclination segments connect pairs of 
the four vertices (Fig. Experimentally, rewiring is 

achieved by local laser melting of a nematic [T3] . 

The three- fold symmetry of —1/2 disclinations (Fig. 2| 
entering the tetrahedron through the vertices coincides 
with the Cs tetrahedral symmetry axes. The director 
field inside the tetrahedron has intrinsic dihedral sym- 
metry (D2d), ensured by the relative positioning of the 
disclinations. Due to this symmetry and the profile of the 
disclinations, the director field stands perpendicularly to 
all the faces of the tetrahedron and makes hyperbolic 
turns at all the edges, thus completing the full tetra- 
hedral symmetry of the director field on the surface of 
the tetrahedron. Consequently, rotations from the tetra- 
hedral symmetry group preserve the continuity of the 
disclination lines and the surrounding director field and 
therefore always generate physically possible structures. 
As the disclination segments inside the tetrahedron have 
lower symmetry (£>2d) than the field on its surface, rota- 
tions around a chosen Cs symmetry axis generate 3 dis- 
tinct configurations of disclinations, depicted in Fig. [T]i. 

The real director field deviates from perfect tetrahe- 
dral symmetry in order to accommodate the proximity 
of the particles and to minimize the free energy. How- 
ever, it only differs from the idealization by a continuous 
transformation, so the topological invariants are not af- 
fected. In further derivations, we assume this symmetry 
to be exact. 

The director field surrounding a disclination may ro- 
tate around the disclination line tangent. In a closed 
loop, rotations are restricted by the fact that the direc- 
tor field must be continuous. The loop, together with the 
orientation of its cross-section, can be described mathe- 
matically by a ribbon (Fig. A ribbon can be assigned 
a self-linking number, SI, a topological invariant that la- 
bels how many times it turns around its tangent in the 
course of one loop. Because of the three-fold symmetry of 
disclinations, the self-linking number is not restricted to 




FIG. 2. To describe a disclination line, we introduce a ribbon, 
denned by an axis curve and a secondary curve that follows 
the orientation of the field cross section of the disclination. 
The ribbon may reconnect with itself with an offset angle of 
±120° while keeping the director field continuous, due to the 
symmetry of its cross section. 



integers, but can assume any third-integer value (Fig. [2]). 
Using of Cdlugdreanu theorem, we can decompose the 
self-linking number into writhe and twist [20j [21] , 



SI = Wr + Tw. 



(1) 



Writhe depends on how the loop changes direction in 
space, while twist contains information about the local 
torsion of the ribbon around its axis. Consider the theta 
structure (Fig. [T^i). Up to an arbitrary homotopic trans- 
formation, the structure is completely symmetric and 
both loops are planar. Both twist and writhe therefore 
equal zero, which can be verified using the corresponding 
Gauss integral definitions [20 . Tetrahedral rotation of 
a portion of the ribbon does not change the twist, as it 
is defined as an integral of local twist density, which is 
preserved by rigid transformations. As we have shown 
that the theta structure has zero twist, the same holds 
for all entangled dimer structures. It follows from Eq. ([!]) 
that in this idealization the self-linking number equals the 
writhe. The physics of nematic liquid crystals is hidden 
in the transformation rules for rewiring and is not in- 
volved in the computation of writhe, which only depends 
on the disclination loop geometry. In practice, disclina- 
tion loops do not necessarily have zero twist, but this 
can always be changed by continuous transformations 
that preserve the topology of the structure. The twist 
and writhe convert into each other under such transfor- 
mations, but their sum remains equal to the self-linking 
number, calculated in our idealized case. 

Since the writhe only depends on the axis curve of a 
ribbon, ordinary loops can be used instead of ribbons. 
We use the tantrix representation: the loop is mapped 
to the unit sphere of tangents. Fuller's formula [16] ex- 
presses writhe in terms of the spherical area A enclosed 
by the tangent indicatrix (tantrix) loop. The writhe 
given by this formula has modulo 2 ambiguity because 
full 47r wraps do not change the tantrix loop, 



Wr 



A 
2^ 



1 mod 2. 



(2) 
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FIG. 3. (a) Tetrahedral rotation changes the positions of two disclination line segments, which also changes the curve, traced 
by the tangent on the unit sphere. The spherical area traced by this curve changes by ±4tv/6 (shown in red), which is directly 
related to the ±2/3 change in writhe and consequently, the self-linking number, (b) Schematic depiction of dimer structures 
and transformations between them. A tetrahedral rotation changes the writhe by ±2/3 if parametrizations of initial and final 
structure are the same. By varying the parametrization, we can calculate writhes of all different dimer structures. Depicted 
here are the theta, chiral omega and figure-eight structures and the pair of Saturn rings, with their respective self-linking 
numbers. Structures that consist of two loops are shown with both possible choices of parametrization. Orientations of loop 
parametrizations are indicated by arrows. 

preserving rewiring changes the total writhe by ±2/3 and 
changes the number of loops, n, by one [22 , we can write 
a conservation law 

q n n 

-(^Sl(A i ) + 2^Lk(A i ,A j ))+n = q mod 2. (4) 

i i>j 

Linking numbers are integers with ambiguously defined 
sign and the number of loops may either increase or de- 
crease by one, hence the modulo 2. This relation is a 
generalized conservation of topological charge g[E] and 
reflects the fact that, due to the presence of line defects, 
only the even/odd parity of q is conserved [I2j[23]. The 
interpretation of q as the topological charge can be jus- 
tified with an example. Consider q homeotropic spheri- 
cal particles in a planar nematic cell, each with its own 
Saturn ring loop [2]. Such a system satisfies the above 
equation as it contains n = q unlinked loops with SI = 0. 
Entangled structures can be reached by applying succes- 
sive tetrahedral rotations, which preserve both the left 
side of Eq. Q and the topological charge. 

The derived formalism can be demonstrated using our 
dimer structures. There are two rewiring sites situated 
symmetrically between the colloidal particles. The 3 • 
3 = 9 possible structures consist of one theta structure, 
two equivalent structures with disjoint Saturn rings, two 
equivalent pairs of chiral omega structures and a pair of 
chiral figure-eight structures. The theta structure has 
zero self-linking number and the others are reached by 
successive tetrahedral rotations. The structures with one 
loop have self-linking number ±2/3 while the structures 



The tetrahedral rotations change which pairs of the ver- 
tices are connected by segments of the disclination loop. 
Each loop segment is planar and maps to a great circle 
arc on the tantrix sphere (Fig. [3^l) . The endpoints of 
these arcs are the tangents through the vertices of the 
tetrahedron, which form a square on the tantrix sphere. 
Rewiring induced by the tetrahedral rotation changes 
the enclosed area by ^47r/6, which by Fuller's formula 
([2| corresponds to a ±2/3 change in writhe. This is 
consistent with the restriction of the self-linking num- 
ber to thirds, imposed by the symmetry of —1/2 discli- 
nations (Fig. [2J. It can be shown that the modulo 2 
ambiguity in Eq. ([2| can be dropped in our case (see 
appendix). As the tangents flip sign if the parametriza- 
tion of a curve is reversed, the ±2/3 change in writhe is 
correct only for tetrahedral rotations that preserve con- 
tinuous parametrization of the loops. If this is not the 
case, the change in writhe can be calculated by finding a 
succession of multiple rewirings that result in the same 
structure (Fig. |3Jd). 

We can generalize the notion of writhe for a union of 
two or more loops Ai that may or may not be linked. 
The Gauss integral that defines writhe decomposes into 
writhes of individual loops, which equal self-linking num- 
bers, Sl(Ai), and linking numbers Lk(Ai,Aj) between 
pairs of loops (see appendix) 

Wr(A U-.-UA n ) = J2 Sl(Ai) ± 2 ^ Lk(Ai, A 3 ). (3) 

i i>j 

Combining this with the fact that every parametrization- 
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with two loops have both self-linking numbers equal to 
zero, which agrees with the conservation law Q. The 
results are shown in Fig. [3J3- 

The derived conservation law Q holds for a set of mul- 
tiple linked loops as a whole. Individual constituent loops 
have a self-linking number of the form p/3, where p is 
even if an even number of disclinations pass through the 
loop and odd in the converse case. This can be shown 
by choosing an idealized model of a director field repre- 
senting such loop and determining mathematically which 
element of the fundamental group the given director field 
represents [I2j [23 . The calculation is carried out in ap- 
pendix. 

We have shown that any rewiring of two orthogonally 
crossing —1/2 disclinations is possible, as the topologi- 
cal requirements are satisfied entirely by the changes in 
the self-linking and linking numbers caused by the ap- 
plication of the tetrahedral rotations. In contrast, +1/2 
disclinations cannot form rich entangled structures, as 
they only allow integer self-linking numbers. Under the 
restriction that only —1/2 disclinations are present, the 
self-linking numbers of the loops are topological invari- 
ants, coupled with surrounding topological charges by a 
conservation law. In confined, chiral and field-affected 
environments, however, the type of the disclination pro- 
file may vary between ±1/2 and twist disclinations [24] . 
Our findings do not apply directly to such cases, as the 
self-linking number is ill-defined if the disclination cross- 
section does not have constant symmetry. 

Our work introduces two important advances in the 
theoretical understanding of nematic braids. By com- 
bining the formalism of differential geometry with the 
characteristics of nematic defects, we are able to show 
that the self-linking number is a topological invariant 
of — 1 /2 disclination loops that successfully differentiates 
between the loops and ensures the conservation of topo- 
logical charge. On the other hand, our explanation of 
local rewiring by tetrahedral rotations gives a qualita- 
tive three-dimensional image of disclinations and resolves 
the behavior of director in complex disclination loop net- 
works seen in experiments and simulations [8] [TDJ [13]. 
The richness of entangled structures increases with the 
number of available rewiring sites, which are more abun- 
dant in chiral systems [4] [13] . The rewiring rules classify 
the set of possible braids and allow a transparent design 
of new structures and guidance of their experimental real- 
ization. Provided symmetry- driven rewiring rules similar 
to our tetrahedral rotations exist, our formalism can be 
extended to systems with different line defects, crossing 
geometries [25], or any system with a well-defined self- 
linking number (e.g. loop DNA [TT} [26]). 

We thank T. Lubensky and R. Kamien for helpful dis- 
cussions regarding the topological interpretation of the 
results. We acknowledge support from the Slovenian Re- 
search Agency (research program P 1-0099 and project 
Jl-2335) and the NAMASTE Centre of Excellence. 



I. APPENDIX 

In the following sections, we supplement the main ar- 
ticle with a review of the basic knowledge of differential 
geometry of loops and ribbons needed to describe general 
structures consisting of nematic disclinations. All curves 
are assumed to be regularly parametrized, so that their 
tangents are uniquely defined. 



Appendix A: Gauss map 

Consider two nonintersecting curve segments A and B 
with parameters s and s' . Each pair of points from the 
two segments define a chord. For a given parametriza- 
tion, we can define a Gauss map that assigns a direction 
of the chord to each pair of points (Fig. [4^). This maps 
from the set of chords (chord manifold) to the unit sphere 
[20] [27] 



(8,8') 



r A (s) - r B (s') 
M*)-r B (s')l 



eS 2 



(Al) 



A rectangle in parameter space is mapped to a patch 
on the unit sphere (Fig. |4|. The area of the patch is 
obtained by integrating the Jacobian of the Gauss map, 
which yields a Gauss integral. 

9{A,B) = ±- [ f t A (s)xt B (s'). r f ] ~ rg( f } dsds' 
4ttJJ \r A (s) - r B (s')| 3 

(A2) 

From above definition, it is evident that the Gauss in- 
tegral is commutative and bilinear under union of curve 
segments (Fig. [4^): 

g(A,BuC)=G(A,B) + g(A,C) . (A3) 

If the curves A and B are closed (we will refer to closed 
curves as loops), the parameters are cyclic and therefore 
the chord manifold has no boundary (it has topology of a 
torus). The patch on the unit sphere also has no bound- 
ary, so its area can only be an integer multiple of 47r, de- 
pending on how many times the map wraps the sphere. 
The Gauss integral of two loops defines the linking num- 
ber Lk(A,B) = Q(A,B), which is the number of times 
the loops pass through each other [20] [21] [28] . 

Instead of a simple loop, consider a ribbon, an infinites- 
imally narrow two-dimensional strip, seamlessly closed 
into a loop. A ribbon can be specified by its two bound- 
ary loops: an axis loop and a secondary loop that differs 
from the axis loop by a infinitesimal displacement [27] . 
A ribbon can twist around its tangent, which is suitable 
for description of —1/2 nematic disclination lines. The 
linking number of the axis loop and the secondary loop 
labels the ribbons according to the number of twists in- 
corporated in the loop and is invariant under continuous 
transformations. To differentiate the linking number of a 
ribbon from the linking number of two arbitrary curves, 
we will label it with symbol SI and call it a self-linking 
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FIG. 4. (a) Two curve segments and a few chords spanning between them. Their directions are marked with arrows, (b) For 
two curves in space, a chord between two points is defined for every pair of parameters (s, s'). The chord directions continuously 
vary with the parameters, (c) Directions of chords lie on a patch on the unit sphere. Gauss integral calculates the area of this 
patch. One curve is represented as a union of two segments (red and blue), which manifests as segmentation of the patch on 
the unit sphere. Gauss integral is bilinear because areas of these patches are additive. 



number (this term is used in the field of theoretical ge- 
ometry in a narrower context [29] ). Even though both 
axis loop and secondary loop of the ribbon have to be 
continous for the self-linking number to be defined, we 
can assign fractional self-linking numbers to a disclina- 
tion loop, if the ribbon representing the disclination line 
runs around the disclination loop multiple times. In the 
case of nematics, the three- fold symmetry of —1/2 discli- 
nations restricts the linking number to third-integer val- 
ues. 

The Calugareanu theorem decomposes the Self- 
linking number into twist that can be seen as contribution 
of locally distributed torsion and writhe that measures 
nonplanarity of the curve. The writhe is a Gauss integral 
of the ribbon's axis loop with itself, Wr(A) = Q(A,A). 
The chord manifold in this case is not a torus but an 
annulus, which has two boundaries, so unlike the linking 
number, the writhe can assume any value [27]. The twist, 
on the other hand, is defined as a single integral along 
the loop and the integrand can be interpreted as a local 
twist density [20 . 

TW = 2^ / " X 9sU ^ ds ( A4 ) 

u is a normalized vector, perpendicular to the axis loop, 
that defines the orientation of the ribbon's cross section 
[20] . Writhe and twist are not topological invariants, 
but the self-linking number is, which means that homo- 
topic transformations of the ribbon only convert between 
writhe and twist, preserving their sum. 



Appendix B: Writhe and the tetrahedral rotations 

We have shown that the tetrahedral rotation changes 
the writhe by ±2/3. We used Fuller's formula that 
does not take into account the entire Gauss integral, but 
only maps the boundaries of the chord manifold. This 
technique misses full 4tt wraps of the sphere, so the writhe 
is undetermined up to an integer multiple of 2. With 
a reasoning described below, we show that the change 
of writhe by the tetrahedral rotations is exactly ±2/3, 
without an undefined additive offset. 

We can represent the loop as a union of two loop seg- 
ments inside the tetrahedron A and everything else X 
(Fig. [5^l). Using bilinearity (Eq. |A3[ ) we can expand the 
result to Gauss integrals between combinations of seg- 
ments. 

Wr(A UX)= Q{A, A) ± Q(X, X) ± 2G(A, X) (Bl) 

The tetrahedral rotation rigidly rotates the segments A 
and preserves X, so the only affected term is the last one. 
This term is a mapping from a set of chords spanning be- 
tween A and X to the space of their directions. In our 
case, the tetrahedron also contains nematic director field, 
which means no other disclinations may be present inside 
it. If the rotation is performed continuously, no disclina- 
tion lines are crossed by the rotating segments in the 
process. We can therefore make a nomotopic transfor- 
mation that stretches every part of the loop outside the 
tetrahedron towards the infinity (Fig^). The chords ex- 
tending from the loop segments A to the points at infinity 
X do not change direction under finite movements of the 
loop segments A. The only chords that change are those 
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FIG. 5. (a) A tetrahedral rotation transforms segments A 
into A\ which only influences the chords between these seg- 
ments and the rest of the loop X. (b) We can perform a homo- 
topic transformation that pushes everything but the rewiring 
site to the infinity. In this limit, rewiring for all the curves 
looks exactly the same and the chords connected to the infin- 
ity do not change direction with rewiring. Note that on this 
sketch, the chords differ because the curve is still at a finite 
distance from the rewiring site. 

extending from the loop segments to the straight paths 
connecting the tetrahedral vertices to the infinity. These 
paths are independent from the rest of the loop (X), so 
the change in writhe for tetrahedral rotation is univer- 
sal. The largest angle for which the chords can change 
their direction is 7r/3, which happens for the limiting case 
at the endpoints of the disclination segments. Together 
with the fact that most chords do not change at all, this 
is not enough to change the area of the patch on the 
unit sphere for 47r. Our result that rewiring changes the 
writhe for ±2/3 is therefore correct, without the modulo 
2 ambiguity suggested by Fuller's formula. 



Appendix C: Homotopy of disclination loops with 
nonzero self-linking number 

To investigate the topological properties of loops with 
nonzero self-linking number, we have the liberty to choose 
any configuration with desired self-linking number, as no- 
motopic transformations do not influence the result. We 



choose a planar ring-shaped disclination loop with radius 
1 and parametrize the surrounding space in toroidal co- 
ordinates with the main angle ip and inner angle 0. The 
basis vectors of a crossection perpendicular to the discli- 
nation line are = cos ifje x + sin^e^ and e z . 

We are interested in disclination lines with winding 
number —1/2, which is the number of times the director 
rotates when we encircle the loop. We start by choosing 
director field for the cross section at ip = 0. 

n(0, 0) = cos(-f )e^ =0 + sin(-f )e z (CI) 

To describe a disclination loop with nonzero self-linking 
number SI, this cross section must rotate rigidly 57-times 
when ij; increases to 2tt. By evaluating these rotations, 
we obtain full specification of the director field around 
the disclination loop. 

n(<9,V0 =cos(|(0-35Z^))e^-sin(|(0-35Z^))e a (C2) 

To describe the topology of the director field on a torus, 
we need two winding numbers [23]. The winding number 
of the small toroidal circle we already know: it equals 
— 1/2. The winding number on the great toroidal circle 
is determined by the mapping n(0, -0), where we fixed the 
value of 6 = 0, and measures the homotopy class of the 
director field encircled by our disclination loop. 

n(0,ip) = (cos(|S7^>) cos^, cos(|S7 , 0) sin^, sin(|S7 , 0)) 

(C3) 

Fundamental group of nematic order parameter only has 
two elements: integer winding numbers correspond to 
a director field, nomotopic to defect-free field and half- 
integer, which corresponds to defect lines. We know 
that in the first case, the vector representation of the 
director field is continuous and in the latter case, it has 
a sign discontinuity between n(0, 0) and n(0,27r). As- 
suming SI = p/3, we get the condition cos(p7r) = ±1. 
Self-linking numbers SI — {...,— 1,0,|,|,...} there- 
fore correspond to director field, homotopic to defect- 
free field, which is achieved if even number of disclina- 
tion lines pass through our loop. Self-linking numbers 
SI = {. . . , —1, — |, ^, 1, . . .} arise in case odd number of 
disclination lines pass through our loop. For two loops, 
their linking number measures how many times they pass 
through each other and therefore chooses between two 
distinct sets of possible self-linking numbers described 
above. 
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